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I .  INTRODUCTION 


If  a  plane  shock  wave  in  a  nonlinear  hyperelastic  material  is  obliquely 
incident  upon  a  plane  material  boundary,  then,  provided  that  neither  the  angle 
of  incidence  nor  the  amplitude  is  too  great,  a  unique  reflection  pattern  will 
be  created  and  will  be  composed  of  shocks  and  simple  waves  all  centered  on  the 
point  of  reflection.  Such  a  pattern  may  be  called  regular  reflection.  Two 
previous  papers  by  the  author  describe  the  basic  situation1'2,  and  a  third3 
clarifies  a  point  in  the  linear  theory  that  simultaneously  guarantees  unique 
solutions  in  the  corresponding  nonlinear  theory.  Other  authors  have  dealt  with 
aspects  of  the  reflection  problem  in  isotropic  solids4'5'6. 

It  happens,  however,  that  as  the  amplitude  of  the  incident  wave  increases 
or  as  the  angle  of  incidence  increases,  the  leading  reflected  or  transmitted 
wave  may  move  so  rapidly  that  it  must  overtake  the  incident  wave  at  the  bound¬ 
ary  no  matter  what  its  direction  of  propagation,  and  thus  the  pattern  of 
regular  reflection  will  be  destroyed.  In  linear  elasticity  the  limiting  angle 
of  incidence  for  regular  reflection  to  occur  is  called  the  angle  of  grazing 
incidence  if,  for  that  angle,  the  incident  wave  is  the  fastest  that  can  occur, 
and  it  is  called  a  critical  angle  if  the  incident  wave  is  one  of  the  slower 
waves.  In  nonlinear  elasticity  the  distinction  seems  artificial  so  the  phrase 
"limiting  angle"  will  be  used  to  cover  all  cases. 

The  essentials  for  shock  reflection  may  be  simply  stated.  For  a  fixed 
angle  of  incidence  the  incident  shock  wave  may  be  specified  completely  by  a 
single  parameter,  thus  fixing  the  deformation  gradient  and  particle  velocity 
immediately  behind  the  wave.  The  problem  now  is  to  fit  reflected  waves  so  as 
to  connect  the  state  just  fixed  with  some  state  at  the  boundary  that  is  com¬ 
patible  with  the  boundary  conditions.  Evidently,  the  situation  is  mathemat¬ 
ically  similar  to  the  Riemann  initial  value  problem7  or  to  the  problem  of  waves 


1T.  W.  Wright ,  Reflection  of  Oblique  Shock  Waves  in  Elastic  Solids ,  Int.  J. 
Solids  Structures ,  7,  161-181  (1971). 

2T.  W.  Wright ,  Uniqueness  of  Shock  Reflection  Patterns  in  Elastic  Solids , 

Arch.  Rat.  Mech.  Anal.,  42,  116-127  (1971). 

3T.  W.  Wright,  A  Note  on  Oblique  Reflections  in  Elastic  Crystals,  Quart.  J. 
Mech.  Appl.  Math.,  29,  16-24  (1976). 

4  (7.  Duvaut,  Phenomenes  De  Reflexion,  Refraction,  Intersection  d  Ondes  Planes 
Uniformes  dans  des  Materiaux  Elastiques  Non  Lineaires,  C.  R.  Acad.  So.  Paris, 
SSrie  A,  264,  883-886  (1967). 

5C.  Duvaut,  Ondes  dans  des  MatSriaux  de  Type  Harmonique.  Reflexion  Oblique 
d  une  Onde  de  Choc  Plane  Longitudinale  sur  une  Paroi  Fixe,  C.  R.  Acad.  Sci. 
Paris,  SSrie  A,  266,  246-249  (1968). 

6S.  R.  Reid,  The  Influence  of  Nonlinearity  Upon  the  Reflection  of  Finite 
Amplitude  Shock  Waves  in  Elastodynamics ,  Quart.  J.  Mech.  Appl.  Math.,  26, 
185-206  (1972). 

7P.  Lax,  Hyperbolic  Systems  of  Conservation  Laws  II,  Comm.  Pure  Appl.  Math., 
10.  537-666  (1957). 


7 


initiated  by  uniform  impulsive  loads  on  the  boundary  of  a  semi-infinite  half 
space8*9.  Although  there  are  differences  in  detail,  most  notably  in  the 
criteria  for  admissibility  of  reflected  waves  and  in  the  condition  that 
reflected  waves  all  have  a  common  reflection  point  rather  than  a  fixed  direc¬ 
tion  of  propagation,  the  methods  for  constructing  a  solution  to  the  reflection 
problem  and  for  establishing  existence  and  uniqueness  for  small  but  finite 
amplitudes  follow  closely  those  used  by  Lax7,  and  are  described  in  detail 
elsewhere2 . 

In  summary,  in  Reference  2  it  was  shown  that  each  reflected  wave  connects 
a  fixed  state  ahead  of  the  wave  with  a  one  parameter  family  of  states  behind 
the  wave.  Furthermore,  parameterization  may  be  so  arranged  that  positive  values 
correspond  to  simple  waves  and  negative  values  correspond  to  shock  waves.  In 
anisotropic  solids  there  are  three  possible  families  of  reflected  waves  so  that 
a  sequence  of  such  waves  connects  the  state  behind  the  incident  shock  with  a 
three  parameter  family  of  states  adjacent  to  the  boundary.  In  general,  there 
will  be  three  independent  boundary  conditions  from  which  to  find  the  parameters 
of  the  reflected  waves  (y^,  y^,  y^)  in  terms  of  the  parameter  of  the  incident 

wave  y  .  Symbolically,  we  have 

Bm(Yl’  Y2*  y3’  Yo)  =  2»  3  C1-1) 

The  situation  for  a  reflection/transmission  problem  is  exactly  analogous,  but 
now  there  will  be  three  transmitted  waves  as  well  as  three  reflected  waves  so 
there  will  be  six  parameters  to  be  found  from  six  boundary  conditions.  In 
special  cases  of  material  symmetry  or  degeneracy  there  may  actually  be  less 
than  three  waves  in  a  reflection  problem  or  less  than  six  waves  in  a  reflection/ 
transmission  problem.  Some  cases  of  this  kind  have  been  discussed  in  Refer¬ 
ence  1 . 

A  simple  application  of  the  implicit  function  theorem  shows  that  (1.1) 
can  be  inverted  to  obtain  three  functions  Ym  =  rm(T0);  m  =  1,  2,  3  provided 

that  yq  is  suitably  small  and  provided  that  the  corresponding  linear  problem 

can  be  solved1*2,  for  the  Jacobian  determinant  of  (1.1)  is  identically  the 
same  as  the  one  that  occurs  in  the  equivalent  linear  problem. 

Here  the  reflection  problem  diverges  sharply  from  the  Riemann  or  the  impul¬ 
sive  loading  problems.  In  the  latter  two  cases,  if  the  material  is  hyperelastic 
with  nonvanishing  wave  speeds,  the  determinant  is  always  nonzero.  However,  in 
reflection  problems  under  the  same  conditions,  it  was  shown  in  Reference  3  that 
the  linear  determinant  is  guaranteed  to  be  nonzero  only  if  the  incident  angle 
is  neither  critical  nor  of  grazing  incidence.  In  these  cases  of  limiting  angle 
no  information  can  be  given  in  general.  Each  case  must  be  examined  separately 
since  examples  show  that  sometimes  the  determinant  vanishes  and  sometimes  not3. 
In  any  event  there  is  always  some  loss  of  differentiability  in  the  mathematical 

®L.  Davison,  Propagation  of  Plane  Waves  of  Finite  Amplitude  in  Elastic  Solids, 

■J.  Meah.  Phys.  Sol.,  14,  249-270  11966). 

9A.  S.  Abou-Sayed  and  R.  J.  Clifton,  Analysis  of  Combined  Pressure-Shear  Waves 
in  an  Elastic, /Vieooplaatic  Material,  J.  Appl.  Mech.,  44,  79-84  (1977). 
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structure  of  the  reflected  waves  as  the  limiting  angle  is  approached.  Nothing 
comparable  occurs  in  either  the  Riemann  or  impulsive  loading  problem. 


This  loss  in  uniformity  of  the  solution  may  be  clearly  seen  in  the  struc¬ 
ture  of  reflected  waves  as  follows.  In  the  construction  of  a  composite  family 
of  simple  waves  and  shock  waves  (see  Reference  2)  the  two  halves  of  the  family 
are  joined  together  at  zero  amplitude  where  the  following  equation  holds. 


4p  V(0)  G*  CO)  |  b (0)  •  t(0)  -  Vh(0)  cos  0(0)} 
=  C3(0)  |r(0)  fi  N(0)}3 


(1.2) 


The  angle  of  a  reflected  shock  is  0,  and  the  prime  represents  differentiation 
with  respect  to  the  parameter  of  the  wave  family.  The  right  hand  side  is  a 
contraction  of  third  order  elasticities  and  is  generally  nonzero.  The  term  in 
brackets  on  the  left  hand  side  has  a  geometric  interpretation.  Terminology 
and  symbols  are  explained  precisely  in  the  sequel,  but  for  now  it  is  enough  to 
comment  that  as  the  limiting  angle  is  approached,  this  left  hand  factor  tends 
to  zero  for  the  leading  reflected  wave  so  that  0'(O)  becomes  infinite. 


In  this  paper  the  reflection  of  weak  but  finite  amplitude  shocks  is 
examined,  with  particular  care  being  taken  to  obtain  solutions  that  remain 
valid  all  the  way  to  the  limiting  angle.  In  the  process,  the  limiting  angle 
itself  as  a  function  of  incident  amplitude  is  obtained.  The  procedure  to 
obtain  these  results  is  one  of  expansions  in  small  parameters  for  a  reflection 
pattern  that  is  close  to  the  pattern  appropriate  for  the  corresponding  linear¬ 
ized  problem.  In  the  next  section  the  structure  of  weak  shocks  and  simple 
waves,  taken  individually  either  as  incident  or  as  reflected  waves,  is  devel¬ 
oped.  In  Section  III  general  procedures  for  solving  reflection  problems  in 
materials  of  any  symmetry  are  laid  out  in  a  step  by  step  fashion.  In  Section 
IV  the  general  procedures  are  used  to  obtain  explicit  solutions  in  isotropic 
materials  for  incident  shock  waves  that  are  either  longitudinal  or  quasi- 
transverse  and  for  several  boundary  conditions  in  each  case. 


i 

i 


l. 


i 


II.  STATEMENT  OF  THE  PROBLEM  AND  WAVE  ANALYSIS 

Suppose  a  plane  shock  wave  with  weak  but  finite  amplitude  is  incident 
upon  a  plane  boundary  in  a  nonlinear,  simple,  hyperelastic  solid,  and  suppose 
that  a  regular  reflection  pattern  of  shocks  and  simple  waves  is  established. 
The  situation  is  shown  schematically  in  Figure  1.  The  deformation  is  given  by 
a  function 


*  =  S(X,  t)  (2.1) 

where  x  is  the  spatial  position  of  a  material  particle,  X  is  the  position  of 
the  same  particle  in  a  fixed  reference  configuration,  and  t  is  time.  Both  x 
and  X  are  referred  to  the  same  fixed  Cartesian  coordinate  system.  The  Piola- 
Kirchhoff  stress  tensor,  T,  is  derivable  from  a  stored  energy  function  W, 
which  in  turn  is  assumed  to  be  a  function  only  of  the  deformation  gradient,  F 
If  the  material  does  not  conduct  heat,  no  conclusions  are  changed  by  letting 
W  depend  on  specific  entropy  as  well  since  entropy  change  is  at  most  third 
order  in  the  amplitude  of  a  weak  shock. 
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Figure  1.  Reflection  Geometry  for  the  General  Case 


W  «  W(F) 

T  =  If  =  TCF)  ;  or  T.a  =  Tia(F.e)  (2.2) 


F  = 


3x 

3X 


or 


Fje 


Hi  = 
3Xe 


j,6 


When  subscript  notation  is  used,  Latin  indices  will  refer  to  spatial  coordi¬ 
nates  and  Greek  indices  will  refer  to  material  coordinates. 


The  relations  to  be  satisfied  across  a  shock  wave  are  as  follows* 0. 


[T]N  +  pv[u]  =  0 

[F]  =  -a  8  N  (2.3) 

[u]  =  V  a 


The  square  brackets  indicate  the  jump  of  a  quantity  across  a  shock;  for  exam- 
ple,  [T]  =  T  -  T  where  the  (+)  and  (-)  signs  refer  to  the  front  and  the  rear 
of  the  shock,  respectively.  The  unit  normal  and  normal  speed  of  the  shock  in 
the  reference  configuration  are  N  and  V.  The  density  in  the  reference  config¬ 
uration  is  p,  the  particle  velocity  is  u  *  ff-  ,  and  a,  defined  by  (2.3).,  is 
the  amplitude.  31  3 

In  a  simple  wave,  all  field  quantities  are  functions  of  a  single  parameter, 
say  y,  and  all  quantities  are  constant  on  propagating  planes,  called  wavelets, 
which  are  also  parameterized  by  Y .  Throughout  a  simple  wave  the  following 
equations  must  be  satisfied11. 


(9  -  PV21)W  =  9 

E'  =  -  m  8  N  (2.4) 

y '  =  Vm 

Here  Q  is  the  acoustic  tensor  with  components  Q. .  =  C.  .  N  N  where 
2  ^3  J  p  ot  p 

j  =  3  W/aF^SFjg.  The  prime  indicates  differentiation  with  respect  to 

Y,  N  and  V  are  the  unit  normal  and  normal  speed  of  a  wavelet  in  the  reference 
configuration,  and  m  is  a  unit  vector.  As  in  the  shock  case  p  is  the  reference 

2 

density.  Naturally,  it  is  required  that  pV  be  a  proper  number  and  m  the 
corresponding  proper  vector  of  Q. 


*0(7.  Trueedell  and  R.  Toupin ,  The  Classical  Field  Theories ,  Flugge's  Handbuah 
der  Physik  III/l,  Springer ,  Berlin- Gottingen-Heidelberg  (1960). 

u£\  Varley ,  Simple  Waves  in  General  Elastic  Materials,  Arch.  Rat.  Mech.  Anal., 
20,  309-328  (1965). 


11 


The  response  to  an  incident  shock  of  weak  but  finite  amplitude  will  be 
near  to  the  response  for  the  corresponding  linearized  problem.  Accordingly, 
all  calculations  may  be  based  on  variables  that  measure  small  variations  from 
the  pattern  and  response  from  infinitesimal  elasticity.  In  Figure  2  the 
reflection  pattern  for  the  linear  problem  is  shown  by  solid  lines,  the  actual 
pattern  is  indicated  by  the  dashed  lines,  and  the  angles  o^,  a2,  c*3  are  small. 

Since  all  waves  are  weak,  the  general  constitutive  equation  (2.2).  need  only 
be  approximated  to  second  order  as  follows. 

T  =  C2  (F-l)  +  Jj  C3  (F-l) 2  (2 . 5) 


C«  and  C,  are  tensors  of  second  and  third  order  elastic  constants, 

2  2  aw 

C~  =  a  W/aF  ,  C.  =  *—=•  where  both  tensors  now  refer  to  an  unstressed  reference 

3f 

state  wherein  F  =  1.  In  component  form,  the  symbol  £L(F-i)  is  the  contraction 

2  ^ 

C.  -0(F.0-<S.  )  and  similarly  C-(F-l)  is  the  contraction 

1CXJP  JPJP  '  j  ~  ~ 

^iajgkY^jfi'^jfp  ^ky^ky^ '  The  more  compact  notation  will  usually  be  used 
hereafter. 


Three  types  of  wave  must  now  be  examined;  the  incident  shock,  a  reflected 
shock,  and  a  reflected  simple  wave. 

(i)  Incident  Shock. 


This  wave  propagates  at  a  fixed  angle  of  incidence  0q,  but  its  speed 


depends  on  the  incident  amplitude  and  will  be  slightly  greater  than  the  corre¬ 
sponding  speed  in  linear  elasticity.  For  the  fixed  direction  of  propagation 
there  exists  in  the  linear  theory  a  tensor  Q  with  unit  proper  vectors  p,  q. 


y  and  corresponding  proper  numbers  pV 


2 

P’ 


We  have 


(Q  -  pVp  J)p  =  0,  etc.  (2.6) 

where  Q  is  computed  as  shown  following  (2.4)  but  using  It  may  be  assumed 

that  the  linear  wave  speeds  are  ordered 


V  >  V  >  V 
p  q  r 


so  that  they  may  be  referred  to  as  fast,  intermediate,  or  slow.  Since  Q  is 
symmetric,  p,  g,  r  are  orthogonal,  and  therefore  the  amplitude  of  the  incident 
shock  has  the  representation 


a  =  cp  +  ng  +  sr  (2.7) 

where  C»  n»  C  are  small  quantities.  The  shock  speed  may  be  written 

VQ  =  V (1+v)  (2.8) 
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Figure  2.  Weak  Shock  Pattern  as  Perturbation  from  the  Linear  Pattern 


where  V  is  one  of  the  speeds  from  (2.6)  and  v  is  a  small  correction.  Finally 
it  may  be  assumed  that  the  deformation  gradient  ahead  of  the  incident  shock, 

F+,  differs  little  from  the  identity  tensor  1,  that  is,  all  components  of  £+-l 
are  small.  The  relationships  among  the  various  small  quantities  are  found  as' 
follows.  With  (2.3)2  an<*  ^*^3  substituted  in  (2.3)^  we  have 

|t(F+  -  a  a  N)  -  T  (F*)}  N  ♦  pV2a  -  0  (2.9) 

Now  put  (2.5),  (2.7)  and  (2.8)  in  (2.9)  with  N  «  Nq  and  resolve  the  resulting 
equation  along  p,  9,  and  r.  To  first  order  the  result  is 

P(V2  -  v2K  ■  0 

p(V2  -  V2)n  =  0  (2.10) 

p(V2  -  V2)5  =  0 


where  use  has  been  made  of  (2.6).  These  three  equations  require  that  V  be 
equal  to  one  of  the  linear  elastic  wave  speeds  (as  already  assumed),  say  V  =  V  . 
Then  the  amplitude  £  is  an  arbitrary  quantity  of  first  order,  but  n  and  C  are  ^ 
zero  to  first  order  provided  V  i  V  and  V  t  V  .  To  second  order  the  compo¬ 
nents  along  p,  g,  and  r  yield  P  ^  P  r 


2pV2  v  -  C3(p  fi  N)2(F\l) 

♦  h  c3(p  a  N)3s  =  0 
c(Vp  -  v2)n  -  c3(p  a  N)  (9  a  N)(F+-in 


(2.11) 


♦  h  c3(p  a  n) 2 (g  a  n)s2  *  0 


There  is  a  third  equation,  similar  to  (2 . 1 1) 2  but  with  the  permutations  q-*-r, 
q-r  and  n-*c .  Equation  (2.11)1  shows  that  v  is  a  first  order  term,  and  (2 . 1 1 ) 2 
shows  that  n  (and  C)  is  a  second  order  term  provided  lVp_Vq|/Vp  >  0(5) ,  i.e. 
provided  and  V  are  not  too  near  to  each  other.  (See  Appendix  A.) 

(ii)  Reflected  Shock. 


A  wave  of  this  type  propagates  at  an  angle  0  +  a,  but  its  speed 
along  the  boundary  is  determined  by  the  incident  wave  at  the  point  of  reflec¬ 
tion.  Refer  to  Figure  3.  The  amplitude  of  the  reflected  wave  is  again  given 
by  (2.7),  repeated  here 


a  =  5p  +  ng  +  ;r  (2.7) 

but  now  p,  g,  r  are  proper  vectors  of  g(N°).  That  is,  they  are  proper  vectors 
of  the  acoustic  tensor  C.  .oN0  N°,  where  N°  is  the  normal  to  the  reflected  wave 

lajp  a  p 

in  the  linear  problem.  Furthermore,  the  speed  is  given  by 
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Figure  3.  Reflected  Shock  Wave 


V  ■  sin (0+o) 


VQ(l+v) 

“slnTT  sin(©+o) 


(2.12) 


and  the  normal  to  the  shock  is  given  by 

8  ■  N°  cos  a  +  t°  sin  o  (2.13) 

where  N°,  t°  are  the  normal  and  the  tangent  to  the  reflected  shock  in  the 
linear  problem. 


Previously,  (2.9)  was  used  to  determine  the  relationships  among  v  and 
other  small  quantities.  Here  we  are  interested  in  the  relationships  among  a 
and  other  small  quantities.  Note  first,  however,  that  the  deformation  gradient 
ahead  of  the  wave  is  determined  by  that  ahead  of  the  incident  shock  plus  the 
sum  of  the  variations  across  all  preceding  waves.  For  example,  if  we  consider 
the  first  reflected  wave,  the  deformation  gradient  ahead  of  it  is 


F+  - 


*0**0 


where  the  subscript  zero  refers  to  the  incident  shock.  It  turns  out,  as  in 
(2.11),  that  the  deformation  gradient  is  only  required  to  first  order  so  that 
only  first  order  terms  need  be  retained  in  aQ.  It  is  convenient  to  write 

simply  F+  in  this  section,  meaning  that  the  sum  of  first  order  variations  are 
to  be  included.  Now  with  (2.5),  (2.7),  (2.12)  and  (2.13)  in  (2.9)  we  arrive 
at  the  desired  equation,  which  may  be  resolved  along  the  proper  vectors  of 

Q(N°) .  The  first  order  results  are  as  follows. 

P  € 


2  -V,2 
P  ho 


sin26 


p  n  V 


sin2® 


P  C 


(v 

(vr2  -  v  2  si„2e) 


(2.14) 


where  V. 


V  /sin  0  . 
o  o 


The  interpretation  of  these  equations  is  similar  to 


that  of  (2.10).  Only  one  of  the  linear  wave  speeds,  say  V  ,  can  equal  sin 

“  o 

In  that  case  €  is  an  arbitrary  quantity  of  first  order,  and  n  and  t  are  of 
higher  order  provided  Vp  /  V^,  Vp  t  Vf. 

To  second  order,  the  components  along  p,  q,  and  r  are  as  follows. 
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2pvp  | b*t  -  vh  cos  e|  o  -  h  c3  (p  a  n)3? 

-2pVp2  v  ♦  C3  (p  fi  N)2CF+-1)  =  0  (2  , 

p  (Vp2  '  Vq2)  n  ‘  -2  *  »)(9  fl  t)  +  (P  8  t)(g  fi  N)}<*£ 

-C3(p  fi  N)(g  fl  tJ)(£+-iK  +  >S  C3(p  fi  N)2(g  fi  N)52  =  0 

There  is  again  a  third  equation,  similar  to  (2 . 15) 2  but  with  the  permutations 

q-*-r,  q+r,  ry»-c  as  before.  The  term  in  braces  in  (2.15)^  may  be  interpreted 

graphically  as  shown  in  Figure  4.  The  vector  b  is  the  radius  of  the  wave 
surface  at  the  point  of  tangency  of  the  reflected  wave.  It  may  be  computed 
from  the  formula 


pV  b  =  C .  . „  p . p . N„ 
v  p  a  iajB  FiF;j  (3 

The  projection  of  b  in  the  X-Y  plane  is  shown  in  the  figure. 

Equations  (2.15)  will  fail  if  either  |V  -V  j/V  or  (b-t  -  cos  0)/V 

p  q  p  0  p 

is  small.  The  first  case  is  treated  in  Appendix  A.  The  second  case  occurs  as 
the  incident  shock  approaches  a  limiting  angle  and  must  be  treated  separately. 
It  is  reasonable  to  assume  that  the  two  cases  do  not  occur  simultaneously  since 
the  fastest  acoustic  speed  is  an  isolated  speed  for  most  solids  and  limiting 
angles  involve  only  the  fastest  reflected  wave.  Thus,  if  nearly  equal  speeds 
occur,  they  will  be  less  than  the  fastest  speed  for  that  direction  and  cannot 
occur  at  a  limiting  angle. 

In  the  case  of  a  limiting  angle,  (2.15).  must  be  replaced  by  an  equation 
2  1 

with  terms  up  to  a  since  the  coefficient  of  a  will  be  a  first  order  term. 

Since  it  is  assumed  that  Vp  and  are  not  nearly  equal,  equation  (2.15)2  is 

sufficiently  complete  for  n  and  neither  (A.2)2  nor  (A.2)3  are  needed.  From 

(2.15)_  it  can  be  seen  that  n  is  0(a£)  at  least,  and  when  insert  id  into  (A. 2)., 
/  2  1 
will  give  rise  to  terms  in  a  .  Considerable  care  must  be  taken  so  that  all 
terms  that  contribute  equally  to  each  power  of  a  will  be  included.  The  equa¬ 
tion  to  replace  (2.15) j  has  the  form 

Aa2  -  Ba  +  C  =  0  (2.16) 

where  A  includes  only  terms  of  order  zero  and  B  and  C  include  first  order 
terms .  We  have 
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Relationship  Between  Wave  Curve  and  Reflection  Point 


B  =  2pV  (b-t  -  V  cos  0)  -  |  C3(p  8  N) 2 (p  8  t)5 

*  n 


+2  C3(p  8  N)(p  8  t) (F+-l)-4pV  Vh  cos  0  v 

V  0 

*C2  [fp  8  N)  (g  8  t)  +  (p  8  t)  (g  8  N)]  . 


2  C3(p  8  N)(g  8  N)  (F  -1)  3  Cj(p  0  B  N)5| 

‘  2 


(2.17) 


2  2 
p(V  -  V  ) 
p  p  q 


2  2 
P<V„  -  V) 


+c2  [cp  8  N)  (r  8  t)  +  (p  8  t)  (r  8  N)] 


2  C3(p  a  N)(r  8  N)(F  -1)  3  C3(p  8  N)^(r  8  N)£ 

2 


P(vp2  -  vr2) 


p(Vp2  -  vr2) 


C  =  H  C3(p  8  M3  +  2pVp2  v  -  C3(p  8  N)2(F+-1) 

Equation  (2.15)^  is  a  specialization  of  (2.16)  when  o  is  a  first  order  quantity, 
(iii)  Reflected  Simple  Wave. 

The  analysis  for  a  reflected  simple  wave  is  slightly  more  compli¬ 
cated  than  that  for  a  reflected  shock  since  the  differential  equations  (2.4) 
must  be  satisfied  throughout  the  wave  rather  than  the  difference  equations 
(2.3)  across  the  wave.  Again  it  is  convenient  to  base  the  analysis  on  the 
proper  vectors  from  (2.6).  Refer  to  Figure  S.  Within  the  simple  wave  the 
acoustic  tensor  Q  depends  on  a  variable  deformation  gradient  and  a  variable 
normal  since 


e2  -  e2  +  e3(F-i) 

fcl  «  y°  cos  a  +  t°  sin  a 
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(2.18) 


Figure  5.  Reflected  Simple  Wave 


where  both  F  and  a  depend  on  the  wave  parameter  y.  The  unit  vector  f ,  which 
is  a  proper  vector  of  Q,  may  be  written 


?  -  +  ng  ♦  «  (2.19) 

where  p,  g,  r,  as  before,  are  the  proper  vectors  of  Q(tj°) .  It  is  assumed  that 
0  is  close  to  one  of  the  other  vectors,  say  m  ■  p,  so  that 

£  =  y l-n2-£2  *  1.  From  (2.4)2  we  may  write 

E(y)  =  E  - J  m(y)  a  N(y)  dy  (2.20) 

With  only  the  lowest  order  terms  retained,  this  becomes 


Y 

f(Y)  =  E+  -  P  0  £J0£y  -  P  0  t°  j  a  dy 

Y  Y 

-  g  a  N°  |  n  dy  -  r  B  N° J 


S  dy 


(2.21) 


Now  with  (2.18),  (2.19),  and  (2.21)  inserted  in  (2.4)r  and  with  V  =  Vh  sin  (0+a) 

a  vector  equation  relating  a,  q,  £  and  y  is  obtained.  At  lowest  order  we  obtain 
(2. IS)  again.  Since  £  =  0(1),  we  have  V  =  V.  sin  0,  as  required  by  linear 

p  o 

theory,  and  n,  C  are  higher  order  terms.  At  first  order  the  following  equations 
are  obtained. 


2pV  (b-t  -  cos  0)a 
*  o 

-c3(p  a  ff°)3  y-  2pVp2  v 

+c  (p  a  N°)2(F+-1)  *  0 

(2.22) 

P(vp2  '  vq2)n  '  C2  [(p  *  n°)(9  a  t°)  +  (p  a  t°)(9  a  N°)]a 
-c3(p  a  N°)(9  a  n°) (F+-i) 

+C3(p  B  N°) 2 (g  B  N°)y  =  0 


There  is  a  third  equation  for  £ ,  similar  to  (2 . 22) 2 -  Provided  that  the 

coefficients  of  a  and  q  are  not  small,  these  equations,  upon  comparison  with 
(2.15),  show  that  to  second  order  a  simple  wave  may  be  replaced  by  a  shock  of 
equal  amplitude  and  located  so  as  to  bisect  the  wedge  of  the  simple  wave.  To 
see  that  this  is  true,  note  that  the  wedge  is  defined  by  a(y)  as  y  varies  in 
the  range  [0,  Yjfor  some  extreme  value  Ym.  With  y  replaced  by  h  y^,  (2.22)  j 
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is  identical  in  form  to  (2.15)j.  Insertion  of  (2.22)2  into  (2.21)  followed  by 

integration  shows  that  f(y )  is  the  same  as  E‘  calculated  from  (2.3),  with  (2.7) 
and  (2.15)2.  m  1 

If  -  Vh  cos  sj/V p  is  small,  then  equation  (2.22) ^  must  be 

replaced  by  one  with  terms  in  a2,  but  (2 . 22) 2  remains  unchanged.  By  the  same 

process  used  to  obtain  (2.22),  but  with  care  taken  to  retain  all  terms  quadratic 
in  small  quantities,  we  again  obtain  an  equation  of  the  form  of  (2.16)  with 
coefficients  A,  B,  C  as  follows. 


A  -  pV  2  cos2©  -  C.(p  I  t)2 


{C2[(p  8  »)(9  *  t)  ♦  (p  *  t)(g  a  N)]}: 
“  ^ 


P(vp2  -  Vq2) 


{c2  [(p  a  N)(r  fl  t)  +  (pa  t)(r  a  n)]}: 

"  ^  n 


2  2 
p(y  -  v  ) 

^  p  r  J 


b  *  2Pvp^-t  -  vh  cos  ©J-2  c3(p  a  N)2(p  a  t)y 

+2  C3(p  a  N)(p  a  t)(F+-l)-4pV  Vh  cos  0  v 

*  o 

+2  C2  [(p  a  N)(g  a  t)  +  (p  a  t)(g  a  $)] 

5  2 

p(V/  -  V  *) 

r  i 

.|c3(p  a  N)(g  a  N)(F+-1)  -  C3(p  a  N)2(g  a  N)y} 
+  2  C2  [(P  a  N)(r  I  t)  +  (p  8  t)  (r  B  N)] 

j.  J 

PCyp  -  V> 

.|c3(p  a  N) (r  a  N)(F+-1)  -  C3(p  a  N)2(r  B  N)y| 
c  -  c3(p  a  N)\  +  2Pvp2  V  -  c3(p  a  N)2(F+-J) 


(2.23) 


Equations  (2.23)  should  be  compared  with  (2.17),  and  (2.22)2  with  (2.15)2. 

obvious  shortened  notation,  the  equations  for  a  in  the  case  of  a  shock  or  a 
simple  wave  are  as  follows. 


Shock  Wave: 

A  o2  -  a  (B2  -  |  B2£  +  b)  +  j  Cj?  +  c  -  0  (2.24) 

Simple  Wave: 

A  a2  -  a  (Bj  -  B2r  *  b)  .  CjT  *  c  =  0 

If  is  large,  then  the  terms  with  A,  B2  and  b  may  be  dropped,  but  near  a 

limiting  angle  B^  will  be  small  so  that  all  terms  must  then  be  retained. 

Strictly  speaking  near  a  limiting  angle  a  simple  wave  can  no  longer  be  replaced 
by  an  equivalent  shock  for  computational  purposes.  To  do  so  would  require  the 
equivalence  of  (2.21)  (withy  =  ?)  and  (2.3)2,  written  here  to  second  order. 


F"  -  F+  -  p  fl  N£  -  p  fl  t  £a 
-cifiNn-rfiN? 


(2. 25) 


If  it  could  be  shown  that  a  good  second  order  approximation  could  be  obtained 
by  setting 

a (y)  dy  (2.26) 


then  comparison  of  (2.22)2  with  (2 .15) 2  would  complete  the  required  equivalence. 

Even  if  (2.26)  does  not  hold  to  second  order,  for  practical  purposes  it  is 
still  useful  and  will  introduce  only  a  small  error.  To  justify  this  last 
statement,  first  solve  (2.24)2  f  or  Y . 


Y 


A a  -  (B^  +  b)a  +  c 
C1  +  B2a 


(2.27) 


The  curve  for  y  as  a  function  of  a  is  nearly  parabolic  for  small  a-  To  repre¬ 
sent  a  simple  wave  there  must  be  a  branch  that  intersects  the  axis  on  which 
Y  =  0  and  that  has  positive  and  increasing  values  of  y  for  increasing  a.  If 
this  condition  has  been  met,  there  are  two  extreme  cases  to  consider.  First, 
if  the  roots  of  the  numerator  of  (2.27)  are  widely  separated,  then  the  branches 
may  be  approximated  locally  near  y  =  0  by  linear  functions.  This  occurs  only 
when  Bj  is  large  and  is  the  case  already  considered.  In  the  other  extreme 

condition  the  roots  of  the  numerator  of  (2.27)  are  equal,  say  a  =  o,  and  we 
have  the  following  approximate  expressions. 


o  -  a  =  Y* 

®*t  ■ a  £  *  j  {3/2 


On  the  other  hand  from  (2.24)j  in  the  same  circumstances  we  have  approximately 
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Since  2/3  and  1//2  differ  by  less  than  6%,  the  approximation  (2.26)  should 
prove  useful  even  in  extreme  circumstances. 


III.  REFLECTION  CALCULATIONS:  GENERAL  CONSIDERATIONS 


All  the  usual  boundary  conditions  require  the  calculation  of  stresses  or 
velocities  or  both  behind  the  reflection  pattern.  With  reference  to  Figure  2, 
the  deformation  gradient  and  velocity  behind  the  pattern  may  be  written: 


h 


3 

E 

i=0 


A. 

l 


(3.1) 


u„ 


*  E 

i=0 


y 


where  the  subscript  B  indicates  the  region  behind  the  wave  pattern  and  adjacent 
to  the  boundary,  the  subscript  I  indicates  the  region  ahead  of  the  incident 

wave,  and  the  symbol  indicates  the  variation  across  the  ith  wave.  In  the 

treatment  of  individual  waves  in  the  last  section,  it  was  shown  that  to  second 
order  all  waves  may  be  treated  as  if  they  were  shock  waves.  (The  only  possible 
exception  is  the  last  case  treated,  i.e.,  a  simple  wave  near  the  limiting 
angle.)  Accordingly,  the  amplitudes  of  the  waves  may  be  represented  as  follows. 


»o  8  5o  Po  +  no  9o  +  ?o  5o 

§1  =  *  n2  gj  ♦  Tj 

§2  =  h  e2  +  +  C2  -2 

*3  =  h  53  +  n3  93  +  C3C1+a))-3 


(3.2) 


In  (3.2)  the  subscripts  refer  to  the  various  waves  as  shown  in  Figure  2.  The 
polarizations  p  ,  a  ,  j  correspond  to  the  directions  of  waves  in  the  linear 
pattern.  Thus  n  n 
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9n 


?n 

% 


=  p 


(3.2) 


where  <0  i  .  .  =  C.  The  normal  to  the  incident  wave  is  given  as  well  as 

\*n(  ij  ioj8  a  6 

the  primary  amplitude  and  polarization  of  the  incident  wave,  i.e.,  one  of 


nQ  or  is  given. 


The  primary  amplitudes  of  the  reflected  waves  are 
Cj,  n2  and  All  other  terms,  including  $,  i|i  and  u  represent  higher  order 
corrections.  Transmitted  waves  could  be  represented  analogously. 


The  calculation  of  the  reflection/transmission  pattern  and  of  the  wave 
amplitudes  now  proceeds  in  a  stepwise  manner  once  the  boundary  condition,  angle 
of  incidence,  and  the  polarization  and  amplitude  of  the  incident  shock  are  all 
specified. 

Step  (i)  The  angles  of  reflection/transmission  for  the  linear  problem  are 
computed.  This  corresponds  to  the  fact  that  the  ref lection/ transmission  pattern 
for  an  incident  shock  of  vanishing  amplitude  is  the  same  as  the  pattern  in 
linear  elasticity.  The  pattern  is  computed  by  the  standard  methods,  e.g.,  see 
Musgrave12. 

Step  (ii)  Corrections  to  the  speed  of  the  incident  shock  wave  and  to  the 
polarization  may  be  computed  from  (2.11)  or  (Al)  as  appropriate. 

Step  (iii)  The  amplitudes  of  all  waves,  given  by  (3.2),  are  inserted  in 
the  expressions  for  stresses  and  velocities,  and  these  in  turn  are  inserted  in 
the  boundary  conditions.  The  lowest  order  terms,  considered  separately,  give 
exactly  the  same  expressions  as  would  be  obtained  from  linear  elasticity.  This 
may  be  written  as  follows. 


B 


(3.4) 


Here  B  is  the  3x3  matrix  of  boundary  coefficients  from  linear  elasticity, 
which  operates  on  the  column  matrix  of  primary  reflection  amplitudes,  and  M, 
which  is  a  linear  function  of  5  ,  is  a  column  matrix  that  originates  from  the 

lowest  order  terms  of  the  incident  wave.  For  illustrative  purposes,  only 
reflected  waves  have  been  considered  in  (3.4),  and  the  incident  wave  has  been 
assumed  to  have  primary  polarization  p  .  Equation  (3.4)  may  be  solved  provided 


12W.  J.  P.  Musgrave ,  Crystal  Acoustics ,  Holden-Day ,  San  Francisco  (1970). 
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that  det  B  ^  0.  In  Reference  13  it  is  shown  that,  generally,  det  B  ?  0  except 
possibly  at  a  limiting  angle.  For  the  moment  let  us  assume  that  det  B  is  non¬ 
vanishing  for  all  angles.  The  case  when  det  B  =  0  at  the  limiting  angle 
requires  special  treatment  and  is  discussed  separately  in  Step  (vii) . 

Step  (iv)  Angle  corrections  a^,  c^,  and  amplitude  corrections  n^,  f^, 

2 »  C2»  £3>  n3  may  be  computed  from  (2.15)  or  (A. 2)  in  all  cases  that  do  not 

involve  a  limiting  angle.  If  reflection  occurs  near  a  limiting  angle,  correc¬ 
tions  may  be  found  from  (2.16)  with  (2.17)  (or  with  (2.23)  if  greater  accuracy 
is  desired  in  the  case  of  a  simple  wave) . 

Step  (v)  The  results  of  (ii) ,  (iii),  and  (iv)  are  used  to  express  the 
second  order  terms  of  the  boundary  conditions.  These  may  be  written  as  follows 

$ 

B  n2  =  N  ^5o,  q,  5i#  04)  (3.5) 

“  «3 

Here  B  is  the  same  boundary  matrix  for  the  linear  problem  as  before,  <j>,  iJj,  u 
are  correction  factors  for  the  primary  amplitudes,  and  the  column  vector  N  is 
quadratic  in  first  order  terms  and  linear  in  second  order  terms. 

Step  (vi)  If  the  angle  of  incidence  is  near  the  linear  limiting  angle, 
then  (2.16)  must  be  used  to  compute  the  nonlinear  limiting  angle.  For  shock 
waves  a  limiting  angle  occurs  when  a  becomes  complex  with  (2.17)  used  in  (2.16). 
For  simple  waves  a  limiting  angle  occurs  when  o  becomes  complex  with  (2.23) 
used  in  (2.16).  For  this  last  case  it  is  necessary  to  set  y  =  0  in  (2.23) 
since  the  limit  occurs  when  the  leading  edge  of  the  simple  wave  overtakes  the 
incident  wave.  Note  that  only  Step  (iii),  the  linear  solution,  is  required  to 
compute  the  limiting  angle. 

Step  (vii)  If  the  angle  of  incidence  is  near  the  limiting  angle  and 
det  B  =  0  at  the  linear  limiting  angle,  then  it  is  necessary  to  combine  steps 
(ii)  -  (vi)  into  one  step.  The  boundary  equations  have  then  become  essentially 
nonlinear,  and  the  simplest  approximate  solution  is  obtained  from  a  set  of 
quadratic  equations  rather  than  a  sequence  of  linear  equations. 


IV.  EXAMPLES:  REFLECTIONS  IN  ISOTROPIC  MATERIALS 

The  procedure  given  at  the  end  of  the  last  section  is  perfectly  general 
for  anisotropic  materials.  No  doubt  numerical  results  could  be  produced  by 
computer  for  arbitrary  symmetries,  but  only  isotropic  materials  seem  suffici¬ 
ently  simple  that  explicit  formulas  may  be  derived. 

In  Reference  13  calculations  were  made  for  three  different  boundary  con¬ 
ditions  with  an  incident  longitudinal  shock  wave.  The  results  of  those 

13T.  W.  Wright ,  Oblique  Reflections,  in  Propagation  of  Shook  Wavee  in  Solids, 
ed.  E.  Varley,  AMD  Vol.  17,Mf>ME,  New  York  (1976). 
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calculations  are  summarized  below  since  they  are  illustrative  of  the  kinds  of 
behavior  possible.  In  addition,  new  results  for  incident  shear  waves  are 
derived.  All  waves  will  be  approximated  by  shocks  as  described  in  Section  II. 
Notation  for  elastic  moduli  is  explained  in  detail  in  Appendix  B. 

Longitudinal  Waves 


If  the  incident  wave  in  Figure  1  is  a  longitudinal  wave  progressing  into 
an  unstressed  region,  the  linear  wave  speed  is  Cj  =  /c^/p ,  also 

Po  =  Nq,  nQ  =  CQ  =  0,  and  the  first  order  correction  from  (2.11)1  is 


v  = 


(4.1) 


where  C  =  ^  ^  and  c^  =  X  +  2y,  X, p  being  the  usual  Lame  constants.  For 

most  materials  C  is  negative  so  that  compressive  waves  (i.e.,  SQ>  0)  correspond 
to  shock  waves,  as  is  well  known. 


There  are  only  two  reflected  waves  in  (3.2),  so  a^  =  0.  Also  =  C2  =  °> 
and  p.  =  N.,  q.  =  t.  for  i  =  0,  1,  2.  The  reflection  angles  for  the  linear 

-  -1  (C2  \ 

problem  are  given  by  0j  =  it  -  0q  and  ©2  =  it  -0  where  0  *  sin  sin  0Qj  and 

C2  is  the  linear  elastic  speed  for  shear  waves.  Angular  corrections  and  a 2 
are  found  from  (2.15)^ 

-2  cu  cot0o  +  h  CUj  +  £0) 


-  (C-c)  sin  2©o  50  =  0 


9  ^2  cos0 

z  c.  C11  sin© 

1  o 


.  s  44  r  r 

“2  +  **  ^7  c 


■  h  (C-c)  +  Cj2  +  C44J 


(C0  ♦  Ej)  =  0 


(4.2) 


(4.3) 


Amplitude  corrections  and  ^  are  ^ound  from  (2.15)2 

(c12  +  c44)(nj_  -  +  h  (C-c)  sin  40Q5OC1  =  0 

-  (c12  ♦  c44)(C2  +  a2n2)  +  **  [*a  (C-c)  +  c12  +  c44]  n2 

+  J*(C-c)  |cosin  2(©0  +  0) -CjSin  2(0Q  -  §)J  n2  =  0 

In  (4.2)  and  (4.3)  primary  amplitudes  £j,  ri2  are  t0  be  found  from  the  linear 
boundary  value  problem,  and  c  =  C1J2  =  cn  \\  22' 
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Note  that  (4.2)^  clearly  fails  as  0Q  approaches  90°,  the  limiting  angle 

in  this  case,  since  the  coefficient  of  a,  tends  to  zero.  With  6  ■  n/2  -  a  , 

1  0  0 

where  aQ  is  small,  (4 . 2) .  must  be  replaced  by  a  quadratic  equation,  obtained 
from  (2.16). 

C11  °12  +  2  C11  °o“l  '  **  C(5o  +  "  0  C4-4) 

Since  must  be  real,  the  discriminant  of  (4.4)  must  be  positive.  It  is  this 

requirement  that  determines  the  limiting  angle  as  a  function  of  amplitude.  The 
discriminant  vanishes  when 


and  then  *  -  aQ  so  that  ©1  =  y  .  *n  (4.5)  the  ratio  Cj/ is  a  func¬ 

tion  of  aQ  and  comes  from  the  reflection  problem  for  linear  elasticity. 

Equation  (4.5)  was  obtained  with  the  approximation,  made  for  computational 
purposes,  that  amplitudes  may  be  calculated  correctly  to  second  order  from 
formulas  for  shock  waves.  If  the  first  reflected  wave  is  actually  a  simple 
wave,  then  the  regular  reflection  pattern  fails  when  the  leading  edge  of  the 
simple  wave  overrides  the  incident  wave,  not  when  the  approximating  shock  does, 
since  it  is  embedded  within  the  real  simple  wave.  These  remarks  should  be 
made  clear  by  reference  to  Fig.  5.  The  net  result  is  that  for  simple  wave  fans 
it  is  necessary  to  set  £./£  *  0  in  (4.5)  to  obtain  the  correct  limiting  rela¬ 
tionship. 

Formulas  (4.1)  through  (4.4)  are  perfectly  general  for  any  incident  longi¬ 
tudinal  wave  and  any  boundary  condition,  but  (4.5)  gives  the  correct  limiting 
condition  only  if  the  linear  boundary  determinant  does  not  vanish  at  the  limit¬ 
ing  angle. 

The  formulas  for  velocities  and  surface  tractions  behind  the  reflection 
pattern,  correct  to  second  order  in  amplitudes,  are  found  by  combining  (2.5), 
(3.1)  and  (3.2)  for  isotropic  materials. 


Ug  •  Ng  2  ^44  2 

-  *stt0  -  Sin  2e0  ♦  a151  cos  0O  -  n2  sin  ©Q 


d+v)Vh 


2  2 
-rij  sin  ©o  -  (C2  -  a2n2)  ~  cos©  sin©Q 

y»  Jd  2 

JlZfyj —  *  (£0  ♦  C ! )  sin  ©o  -  »sn2  sin  2©  -  +  a^)  sin  2©c 

2-  .  2- 
+  <*2n2  cos  ©  +  C2  sin  © 


r-  »B  •  ™B>  ■  -  «o  *  * 2  c<,!‘9o)-"2  sin  25 

44  \  44 


(nx  ♦  o^j)  sin  2©o  -  ~  (5 2  ♦  a2n2) 

44 


/  2-  2 
-  2(?2  cos  0  +  °2ri2  s*n 


!§) 


2  C  *. 
sin  0  +  - —  cos'-©  ]  +  ?*n 


o  c 


44 


)  •  -4  s  ■&  •■)! 


(4.6) 


+  *0*1 1  b~  (1  +  2  cos2qo)  si"20o  +  I-  C0S  20  o  cos20o  | 
(  44  44  ) 

♦  h  ~n7  cos©  j  sin(en  +  2§)  +  e,sin  (ert  -  2©) 

C44  Z  ° 


l —  JB  •  (™B)  =  (CQ  -  C:)  sin  2©o  -  n2  cos  2©  +  Cnx  ♦  a cos  2©c 


'44 


U2  -  a2n2)  sin  2©  +  h  (tQ2  -  Sj2)  (C-c)  sin 


20, 


-  «w2(e0  *  € 


■>lfe 


+  1  ]  +  Jj  cos  2© 
/  C44 


For  a  clamped  boundary,  the  first  two  of  (4.6)  must  be  identically  zero, 
and  similarly  for  the  last  two  in  the  case  of  a  free  boundary.  The  linear 
boundary  determinants  in  these  two  cases  are  as  follows. 

Clamped: 


det  B  *  cos  (©Q  -  0) 


-  c2/cx  as  ©o  ■+  | 


(4.7), 
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Free: 


det  B 


(—+2  cos2© 
\44 


cos  20+  sin  20^  sin  2© 
o 


12 

— — -  as  © 

C44  C11  0 


£ 

2 


(4.7)  2 


Since  the  determinants  remain  positive  all  the  way  to  the  limiting  angles 
(except  for  a  Poisson's  ratio  of  zero  in  the  second  case),  the  first  five  steps 
at  the  end  of  Section  III  are  all  that  are  needed  to  determine  solutions  to 
second  order.  These  two  cases  were  discussed  further  in  Reference  13. 


The  rigid-lubricated  boundary  by  definition  requires  the  first  and  the 
fourth  of  (4.6)  to  vanish.  The  linear  terms  give  the  solution  =  1, 

^2l^0  =  0*  but  the  boundary  determinant  in  that  case  is  given  by 


det  B  =  sin©  cos© 
o  o 

+  0  as  S 

o  2 


(4.8) 


Since  sin©Q  appears  in  (4.8)  only  because  of  the  normalization  by  Vh  in  (4.6) 

there  is  no  difficulty  at  0q  =  0,  but  as  ©Q  approaches  90*  the  linear  solution 

no  longer  gives  a  good  first  approximation.  It  is  then  necessary  to  use  non¬ 
linear  boundary  conditions  right  from  the  start. 

With  0q  =  j  -  aQ,  where  aQ  is  a  small  angle,  and  with  a  small  angle 
approximation  used  throughout,  it  is  readily  found  from  (4 . 2) 2  and  (4.3)2  that 


«2  -  o  «0) 

^2  =  ®  ^a2n2^ ' 

Furthermore,  from  (4.3)^  and  (4.4)  we  have  to  lowest  order 

nl  =  ai*i 

= -•  *  h2  * 15  ni «• 4 

Boundary  condition  (4.6) ^  reduces  to 

c44 

7^  *2  “  (5o  '  *15  °o  "  “1*1 


(4.9) 


(4.10) 


(4.11) 
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where  (4.10)^  has  been  used.  Finally,  boundary  condition  (4.6)^,  with  the  use 
of  (4.9),  (4.10)  and  (4.11)  becomes 


Ci  r  c  c0  /  c^-ft 

1  ♦  r  *  ♦  — V  1  +  r)  - 0 
5°  L  2cn“o  [  Wj 


(4.12) 


This  equation  has  only  one  acceptable  real  solution. 

Ci  /  ^  1^5 


c* 15  d2-°}' 


11  “0 

and  4  i  D  <  • 


(4.13) 


The  extreme  case  of  D  large  leads  to  the  linear  elastic  result,  but  the  limit 
D  *  4  gives 


&•* 


(4.14) 


and  the  limiting  amplitude-angle  relation 

£o  *  '  2C7c 


“o2 


11 


(4.15) 


The  results  in  (4.14)  and  (4.15)  were  previously  given  by  Reid6,  but  not  (4.13). 
Note  that  (4.15)  is  similar  to  (4.5),  but  much  stronger. 

Equation  (4.13)  is  plotted  in  Fig.  6.  The  limiting  angle/ amplitude  rela¬ 
tions,  given  by  equation  (4.5)  for  a  simple  wave  and  by  (4.15)  in  the  present 
case,  are  plotted  in  Fig.  7  for  elastic  constants  C/c^  =  -8.5.  Fig.  6  shows 

the  rapid  change  in  reflection  amplitude  as  the  limiting  angle  is  approached. 
Fig.  7  shows  how  large  the  limiting  angle  ac  may  be  even  for  relatively  small 

strains.  This  point  is  discussed  further  in  Reference  13. 

Shear  Waves 


If  the  incident  wave  in  Fig.  1  is  a  shear  wave  progressing  into  an 
unstressed  region,  the  linear  elastic  wave  speed  is  =  /c^/p  ,  with  pQ  =  tQ, 

and  first  order  corrections  v  and  from  (2.11)  are  given  by 


V  a  0 

(4.16) 

-(c12  +  c44}  Co  +  f  {C_c  +  2  (c12  +  C44}}  %2  “  0 
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Figure  6.  Reflected  Shock  Amplitude  at  a  Rigid/Lubricated  Boundary  for  an  Incident 
Longitudinal  Shock.  (Symbols  are  explained  in  the  text) 


Again,  there  are  only  two  reflected  waves  with  angles  ©2  "  ”  ®o  aiul  ®l  *  v  '  ® 


where  5 


■  sin"1  (  —  sine  \ 
-1  C1  '  2  ' 


The  linear  critical  angle  occurs  for  0  ■  y  or 


sin  —  .  Angular  corrections  and  02  are  found  from  (2.15)^. 


2  CU  if  HST  “1  -  *  C  'l  -  ‘“C-clsin  W,  -  Vo  ■  0 
1  0 

C°seo  r 

2  C44  slne^  °2  "  *  LC'C  +  2(C12  *  C44\J  51  *  0 
Amplitude  corrections  n ^  and  £2  are  found  from  C2 - 15) 2 * 


(4.17) 


(C12  +  C44)(nl  "  al5P 

+  H  £(C-c)cos  2(eo  +  0)  +  2  (c12  +  c44)]  n05],  =  0 
(c12  +  C44HC2  *  a2n2^  (4-18) 

-  !*  £(C-c)cos  4  0q  +  2  (c j 2  +  c44)]  non2 
+  h  (C-c)sin  2(0  -  0o)^1n2  -  |  [c-c  +  2(c12  +  c44)J  n2 2  *  0 


As  before  ^  and  r)2  are  t0  t>e  found  from  the  linear  boundary  value  problem. 


rhe  equation  for  again  fails  near  the  limiting  angle  since  the 

-1  c2 

coefficient  of  a  in  (4.17).  tends  to  zero.  With  0  *  sin  - -  a  and 

X  X  W  v  ^  U 

®1  +  “l  *  T  +  “1*  w^ere  a0  an<*  are  small  angles,  eqn.  (4.17)j  must  be 

replaced  by  a  quadratic  equation,  obtained  from  (2.16)  as  before.  It  is  more 
convenient  to  work  with  6^  rather  than  a^.  To  be  consistent  the  basis  vectors 

for  the  first  reflected  wave  now  are  *  ex  and  g^  ■  Eqns.  (4.18)  remain 

correct  with  5  *  ^  and  5^  and  interpreted  as  components  along  the  redefined 

Pl  and  gl. 


-  2  .  12 
:11  “1  +  CC"c)  c~  Vl 


-H  C  ( 


l 


) 


(4.19) 
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The  last  term  in  (4.19)  does  not  cone  from  (2.17)^,  but  rather  from  (2.14) ^ 
since  the  use  of  instead  of  aj  leaves  a  second  order  term  left  over  from 


(2.14)^  that  must  be  accounted  for. 


*1  =  - 


Solution  of  (4.19)  gives 
2 


\  U 


(4.20) 


where  the  positive  radical  has  been  chosen  so  that  6^  has  the  correct  value 


for  fixed  a  with  n  =0. 
o  o 

radical  must  be  positive. 


To  ensure  that  is  real,  the  terms  under  the 


As  before,  it  is  this  condition  that  determines  the 
limiting  relationship  between  nQ  and  aQ. 


The  ratio  must  be  obtained 


from  the  linear  problem  (or  nonlinear  problem)  if  denotes  a  shock  wave  or 
be  set  equal  to  zero  if  denotes  a  simple  wave  fan.  It  should  be  noted  that 
in  the  development  given  here  is  required  to  be  positive  since  only  then 
can  C^nQ  be  found  from  the  linear  problem.  If  <  0,  perhaps  an  estimate 
for  a  lower  bound  of  nQ  could  be  found  in  some  cases  by  appealing  directly  to 
the  nonlinear  problem.  Of  the  two  values  of  tiq  that  make  the  radical  be  zero 
only  the  smaller  one  is  physically  reasonable. 


(4.21) 


It  is  worth  considering  (4.20)  and  (4.21)  further  for  clamped  or  free 
boundary  conditions.  In  these  cases  it  is  readily  found  that  the  linear  solu¬ 
tion  gives 


Clamped: 


sin  20. 


1 


n  c. 

cos  (9  -  0  ) 

c2 


sin 


1  C2 


(4.22) 
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Free: 


«1 

”o 


sin  40 

o 


cos 


20  „  +  sin  20„  sin  20 
o  o 


(4.23) 


Since  shocks  are  compressive  in  most  materials,  for  the  present  discussion  it 
will  be  assumed  that  C  <  c  <  0.  For  the  clamped  case  if  nQ  >  0  near  the  limit 

ing  angle,  then  <  0  which  indicates  a  simple  wave.  Then  must  be  set 

equal  to  zero  in  (4.21)  and  nQ  has  the  limits 


0  *  %  < 


C11  C1 
C-c 

C2 


(4.24) 


On  the  other  hand, 

n  has  the  limits 
o 


if  nQ  <  0,  then  Cj  >  0,  which  indicates  a  shock  wave,  and 


2 


(4.25) 


If  c  is  actually  positive,  rather  than  negative  as  assumed,  then  (4.25),  or 
more  properly  (4.20),  seems  to  indicate  that  aQ  may  be  negative.  For  the  free 

boundary  case  similar  considerations  apply.  If  nQ  >  0,  then  >  0,  indicating 

a  shock,  and  nQ  is  limited  by  (4.21).  But  if  nQ  <  0,  then  <  0,  indicating 

a  simple  wave,  and  (4.21)  or  (4.20)  seem  to  indicate  that  aQ  may  be  negative, 

-1  c2 

i.e.,  that  the  limiting  angle  is  greater  than  sin  —  .  Fig.  8  shows  an 

C1 

example  of  the  limiting  relations  between  nQ  and  aQ  for  these  two  boundary 
conditions. 


The  boundary  determinants  for  the  clamped  and  free  boundary  have  the  same 
limits  as  (4.7) ^  and  (4.7) 2  so  that  second  order  solutions  may  be  found  by 

following  the  first  five  steps  at  the  end  of  Section  III. 

The  rigid- lubricated  boundary  again  provides  the  most  interesting  case 
since  the  linear  boundary  determinant  vanishes  as  in  (4.8),  and  again  it  is 
necessary  to  use  nonlinear  boundary  conditions  right  from  the  start.  The 
formulas  for  shear  traction  and  normal  velocity  behind  the  reflection  pattern. 
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Figure  8.  Limiting  Angle  -  Amplitude  Behavior  for  Incident  Shear  Waves,  Points  to  the  Right 
of  the  Limiting  Lines  are  Allowed.  (Symbols  are  explained  in  the  text) 


correct  to  second  order  in  amplitudes  are  found,  as  before,  by  combining  (2.5), 
(3.1)  and  (3.2)  for  isotropic  materials. 


-  (n0  ♦  n2)  cos  2®0  -  (50  -  52)  sin  2e0  + 


V^B  •  uB  -  0 


.2 


(n0  ♦  n2)  sin  en  *  Un  -  £,)  sin  0O  cos  e, 


o  *2' 


(4.26) 


-  +  a2  n2  sin  0Q  cos  0Q 

Since  the  linear  problem  has  the  solution  K^j^\0  -  0.  n2/n0  =  _1  for  every 
incident  angle,  it  seems  reasonable  to  assume  here  that  Cj  =  o(nQ)  and 
n2  =  °(n  )»  i.e.,  lim  €,/nQ  *  0  and  lim  n2/nQ  =  constant.  With  these  estimates 
°  "o'*0  °  no'K) 

of  magnitudes,  it  is  easily  seen  from  (4.16)2>  (4 - 17) 2 ,  (4.18)1  and  (4.18)2» 
respectively,  that 

*o  =  0(no2) 


a2 

n 


°(n0) 


(4.27) 


i  ~  ai  h  *  o(no  5 


e2  =  °(n0  ) 

Since  terms  higher  than  quadratic  have  already  been  eliminated,  (4.26)  may  be 
written 


-(n0  +  n2)  cos  2©o  -  UQ  -  e2)  sin  20q  +  2«1C1  =  C 

-(n0  +  n2)  sin20o  ♦  *s(^0  -  K2)  sin  20q  -  OjCj  *  0 
To  second  order  we  have 

n2  -  -n0  ♦  o(n02) 

®l5l  "  (?o  '  ?25  sin0o  cos9o  +  o(no2) 


(4.28) 


(4.29) 
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Equation  (4.29) ^  shows  that  the  reflected  shear  wave  has  nearly  the  same 

amplitude  as  that  calculated  in  the  linear  problem.  However,  (4.29)_  leads  to 
rather  different  results  for  1 


-12 

where  0  =  sin  —  .  For  purposes  of  graphical  interpretation,  (4.30)  will 

c  Cj 

be  rewritten  as 

-r 

Since  the  positive  square  root  is  always  intended,  has  the  same  sign  as  R, 

which  depends  both  on  the  magnitude  of  0  and  the  sign  of  C-c.  Thus  must  be 
a  root  of  ci 


f(?j)  -  A  -  B  5X2  +  R2 

which  is  plotted  in  Fig.  9  for  A  >  0.  In  this  case  if  R  >  0,  the  reflected 
wave  is  a  shock,  so  the  limiting  angle  occurs  when  the  minimum  in  f(Cj)  just 

touches  the  axis.  The  minimum  lies  at  K,  =  2B/3A  and  the  function  vanishes 

2/3  1 

there  if  B  =  3  ft  AR)  .  Returning  to  the  coefficients. of  (4.30),  the  limit¬ 
ing  conditions  are 


On  the  other  hand,  for  A  >  0  but  R  <  0,  the  reflected  wave  is  a  simple 
wave.  With  *  0  for  the  leading  edge  of  the  wave  in  (4.20),  the  radical 
vanishes  at 
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Normalized  to  A  =  R  =1.  (Symbols  are  explained  in  the  text) 


(4.32) 


a 


o 


"o  *  0  i\) 


The  amplitude,  as  calculated  by  the  shock  wave  approximation  from  (4.30),  is 
given  by 


which  is  smaller  than  (4 - 31) 2  by  a  factor  of  4.21^'5. 

The  argument  for  the  case  A  <  0  follows  the  same  pattern  as  for  A  >  0,  but 
with  replacing  at  every  step,  since  for  A  <  0,  Fig.  9  becomes  its  own 

mirror  image  by  reflection  across  the  axis  =  0. 


I 
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APPENDIX  A 


NEARLY  EQUAL  WAVE  SPEEDS 

If  Vp  =  Vq  in  (2.11)  or  even  if  they  are  near  to  each  other,  i.e. 

I V  -  V  l/V  «  1,  then  (2.11)-  fails.  In  that  case  n  oust  be  assumed  to  be  a 

1  p  q1  q  2+ 

first  order  term.  With  £  and  F  -1  still  arbitrary  first  order  terms,  it  is 
necessary  to  find  v  and  n  from  simultaneous  quadratic  equations. 

2pVp2  v£  -  C3(pfiN)2(F+-l)£  -  C3(p»N)(g«N)(F+-l)n 

+  h  C3(p«N)3  £2  +  C3(paN)2(q8N)£n 


+  %  C3(p«N)(q«N)2  n2  =  0 
p(Vp2  -  Vq2)n  +  2pV2vn  -  C3(g«N)2(F+-l)n 


(A.  1) 


-  C3(pfiN)(q8N)(F+-l)£  +  H  (p«N) 2 (q®N) £  2 
+  C3(piN)(gaN)2£n  +  h  C3(g#N)3n2  =  0 


The  equation  for  £  remains  unaltered.  Equations  (2.11)  are  specializations  of 
(A.l)  with  n  a  second  order  quantity  and  terms  of  third  or  fourth  order  dis¬ 
carded  . 


For  the  case  of  nearly  equal  speeds  in  (2.15)  it  is  necessary  to  assume 
that  n  is  of  the  same  order  as  £,  and  n  and  a  must  be  found  from  the  follow¬ 
ing  quadratic  equations  where  £,  v  and  F+-l  are  taken  to  be  known. 


-  2pV  (b  •  t  -  V.  cose)a£  +  2pV  2  v5 

P  ^  P  Q  P 

-  C2  {  Cp8N)  (g8i)  ♦  (p«t)(qfl«}an 

-  C3(p*N)2(F+-l)5  -  C3(piN)(9fiN)(F+-l)n 

♦  h  C3(pflN)V  +  C3(pflN)2(gflNHn 

+  H  C3(pfiN) (g|W)2n2  =  0 

P(vp2  -  Vq2)n  -  2p(Vqbq  .  t  -  vpvh  cos0)an  +  2PVp2  vn 

-  C2  { (piN) (git)  +  (pat) (qBN)}  a? 

-  C3(p8N)(q8N)(F+-l)5  -  £3 (gBN) 2 (F+-l)n  (A. 2) 

+  h  C3(paN)2(q8N)C2  +  C3(paN)(qBN)2Cn 

♦  is  c3c9aN)3  n2  -  o 
p(Vp2  -  Vr2)C  -  C2  [(pfiNHrBt)  ♦  (gat)(rBN)]aS 

-  c2  [(gflN) (rat)  ♦  (gat) (riN)]  on 

-  C3(paN)(raN)(F+-l)£;  -  C3(gaN)(rBN)n 

+  h  c3(paN)2(raN)c  2  +  c3(paN)(qay)(raN)5n 

+  h  c3(9aN)2(raN)n  2  =  0 

The  vector  bp  is  the  same  as  ]?  in  the  formula  following  (2. IS),  and  bq  is 

computed  from  the  same  formula  with  the  permutations  p-*q,  p+g.  Equations 
(A. 2)  reduce  to  (2.15)  if  n  is  a  second  order  term  and  terms  higher  than 
second  are  dropped. 

In  case  |V  -  V  |/V  «  1  in  (2.22)  then  n  and  £  are  quantities  of  com- 
P  H  i 

parable  magnitude,  and  both  (2.21)  and  (2.22)  must  be  replaced.  Since  m  is  a 
unit  vector,  the  condition  £  =  1  must  be  replaced  by 
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(A.  3) 


which  is  correct  to  second  order.  The  lowest  order  terns  in  (2.20)  now  are 
given  by  F  «  F*  -  pAN  -  qAN  ryy»  and  by  the  sane  process  used  to  obtain 
(2.24)  we  have 

2pVp(b  •  t  -  Vh  cos0)oC  -  2pVp2  VC 
♦  C2  {(pANHgAj)  +  (git) (gAN)}  on 

♦  C3  (pAN)2(F+-l)C  +  g3(pAN)(qAN)(F+-l)n 
-  C3(p«N)3  c\  -  2  C3(p8N)2(g«N)CnY 

-  C3(pAN)(q«N)2  n2y  =  0 

p(Vq2  "  Vp2)n  +  2p(Vq~q  '  5  "  VpVf  COS0)fl,n  *  2pVp2  vn 

♦  C2  {(pAN)(gAt)  ♦  (pAt)  (gAN)f  <*C 

♦  g  (pAN)(gAN)(F+-l)C  +  C,(gAN)2(F+-»n 

(A. 4) 

-  C3(pAN)2(gAN)  C2Y  -  2  g3(pAN) (gAN)2CnY 

-  g3(gAN)3  n2Y  =  0 

p(Vr2  -  Vp2)c  +  C2  {(pAN)(rAt)  +  (pAt) (rAN)f aC 
+  C2  {(qAN) (rAt)  +  (gAt) (rAN)}  an 

♦  C3(pAN)(rAN)(F+-l)5  +  CjCgAN) (rAN) (F+-l)n 

-  C3(pAN)2(rAN)c2Y  -  2  ?3 (pAN) (gAN) (rAN) Cny 

-  C3(qAN)2(rAN)n2Y  -  0 
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The  interpretation  of  fc  and  ^  is  the  save  as  in  (A. 2).  Clearly  equations 
(A. 4)  reduce  to  (2.22)  when  IVp-V^I/V^  is  not  small,  5*1,  and  n  is  small. 

Equations  (A. 3)  and  (A. 4)  are  to  be  solved  for  a,  5,  n,  c  as  functions  of 
Y«  Comparison  of  (A. 4)  with  (A. 2)  shows  that,  for  purposes  of  calculation,  a 
simple  wave  may  again  be  replaced  by  a  shock  wave.  To  see  that  this  is  true, 
first  note  that  5,  n,  C  in  (A. 2)  are  components  of  a,  which  is  not  a  unit 
vector.  To  make  them  comparable  with  the  terms  of  (A. 4)  we  replace  5,  n,  C 
in  (A. 2)  by  5yb»  nYm*  CYm  where  Ym  ■  |g|.  Next  consider  a  power  series  solu¬ 
tion  of  (A. 4),  5  ■  5q  +  C^y  ♦  •••>  n  ■  +  r»0 Y  ♦  •••><*■  oQ  +  a^Y  +  •••.  The 

terms  50,  nQ,  aQ  are  obtained  from  a  quadratic  set  of  equations,  and  n^, 

from  a  linear  set.  These  results,  when  inserted  into  (2.20),  upon  integration 
give  f(Ym)  up  to  terms  of  second  order  in  small  quantities.  Similarly,  a  power 

series  solution  of  the  revised  form  of  (A. 2)  may  be  sought.  The  leading  terms 
are  the  same  as  those  obtained  from  (2 . 20) ,  but  each  of  the  second  terms  is 
equal  to  one-half  the  value  previously  obtained.  Calculation  of  E(Ym)  from 

(2.3)2>  however,  gives  the  same  result  to  second  order  as  that  obtained  from 

the  assumption  of  a  simple  wave.  It  may  be  concluded  for  this  case  as  before, 
that  a  simple  wave  may  be  replaced  by  a  bisecting  shock  wave  for  second  order 
calculations. 
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APPENDIX  B 


ELASTIC  MODULI 

In  the  literature  on  ultrasonic  wave  propagation,  it  is  commonly  assumed 
that  the  strain  energy  has  the  following  expansion  with  respect  to  an  unstressed 
reference  configuration. 14 


4:<J  ■  i 


Cagy6  ''ag  ^yfi  +  6  cagydyv  ^ag  ^yg  nyv  + 


(B.l) 


Where  nag  =  h  (x^  a  x*  g  -  g^^  is  termed  nonlinear  strain,  and  cagY(5  and 

Cagygyv  are  terme<1>  respectively,  second  and  third  order  tensors  of  elastic 

moduli.  Because  of  symmetries,  the  values  of  the  individual  moduli  are  usually 
reported15  as  entries  in  the  totally  symmetric  arrays  c^  and  ca^c,  where 

a,  b,  c  =  1,  2,  •••*,  6.  Each  Latin  subscript  corresponds  to  a  pair  of  Greek 
subscripts;  1  «-*■  11,  2  -*-+■  22,  5  33,  4  23  or  32,  5  «-*■  13  or  31,  6  ■*-*•  12 

or  21.  Thus,  cn  =  cim,  c23  =  c2233>  c4g  =  c2313,  c15g  =  c111313>  etc. 

The  elasticities  for  use  in  this  paper  may  easily  be  calculated  from  [B-l) . 


'iajg  Xi,0  Xj,p(Caogp  +  Coogpyv  'Vv) 


*  6. .  (c  .  n  *■  h  c  „  n  n  \  +  • 
lj  V  agyv  yv  otgyvpo  yv  pa) 


C, 


iajgky  “  xi,o  Xj,p  CaagpyT 


(B.  2) 


+  6..  x.  c  +  6..  x.  c„ 
ik  j,a  ayga  jk  i,o  gyao 


+  g.  .  x.  c  B  +•••. 
ij  x,o  agyo 

If  the  region  ahead  of  the  incident  shock  wave  is  unstressed,  then  in  (B-2) 

x.  *  5.  and  n  „  =  0.  The  elasticities  reduce  to 
i9a  la  ap 


14tf.  Brugger,  Thermodynamic  Definition  of  Higher  Order  Elastic  Coefficients , 
Phye.  Rev.,  133,  A16U-A1612  (1964). 

15J?.  Beckman  and  R.  F.  S.  Hearman,  The  Third  Order  Elastic  Constants,  in 
Landolt-Bdmstein,  Numerical  Data  and  Functional  Relationships  in  Science 
and  Technology ,  New  Series ,  Group  III ,  Vol.  2,  e<f.  K.-H.  Reltwege,  Springer- 
Verlag,  Berlin ,  Heidelberg ,  New  York  (1969). 


^iajB  "  jlo  *jp  cao8p 


CiajBkY  “  aio  *jp  ^kt  caaBpYT 


4  *ik  ^jo  c#y8o  4  3jk  ^io  cByoo 


+  6ij  6ka  cafiYff  * 

The  reduced  notation  may  also  be  used  for  £,  and  Thus,  C...  - 

C11  11  11'  C112  "  C11  11  22'  etc*  “ 

For  the  isotropic  case,  there  are  only  two  independent  second  order 
moduli  and  three  independent  third  order  moduli15  and  then  (B-3)  becomes 

CiajB  =  c12  6iu  6jB  +  c44  (6ij  6a6  4  ®iB  6ja) 

CiajgkY  =  C123  5ia  fijB  5kY  +  C144  (fiia  6jy  5kB  4  5iY  6jB  3ka  *  6iB  5ja  ^kv) 

+  (C12  +  c144)  (5ia  fijk  6  By  +  6ik  5jB  &ay  +  fiij  \y  ^afi) 

(B-4 

♦  h  (c.c-  -  c ...  ♦  2  Cj*}  ( 6 , „  6.  6,  +6.  6.  &  . 

\  155  144  44/  V  iB  jy  ka  iy  ja  kB 


+  6,.  6.  6  +6..  6  +  6..  6,  6„  +  6.„  6..  6 

XJ  ka  By  ij  kB  oy  xk  ja  By  iB  jk  aY 

+  6..  6.  S  .  *  6.  6..  6  \  -  Caa  (6.a  6.  ♦  6.  6.  6,  A 

xk  JY  «B  iy  jk  aB/  44  V  iB  jy  ka  iy  ja  kB/ . 

For  isotropic  materials  we  also  have  the  identities 

c44  =  **  (cll  "  c12)’  C144  =  **  (C1 12  "  C123)'  C155  *  *  (Clll  “  c112)‘ 
In  this  paper  the  third  order  moduli  most  commonly  used  are 

C  “  Clll  "  clll  4  3  C11 


(B-5) 


c  "  C112  “  C112  4  c12 
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LIST  OF  SYMBOLS 


A,  B,  C 
a 

B  ,  B 
in 

V  *q 

?2’  >iaj$ 

-3’  ^iajpky 
C,  c 

Caf$y6’  CaBy6yv 

Cab’  Cabc 

Ci*  c2 
D 


e  ,  e 
-x*  ~y 

E-  FiB 

-B’  Ei 

m 

N,  N° 

O(-).  o(-) 

?•  3'  l 
Q 

Tia 
S»  ?° 

t 

u 

V,  V_ 


Vp*  vq’  vr 


Polynominal  coefficients  in  eqn.  (2.16) 

Shock  amplitude  vector,  defined  by  (2.3)^ 

Boundary  equations,  boundary  matrix 

Propagation  vectors,  defined  by  the  formula  following  (2.15) 

Second  order  elastic  moduli 

Third  order  elastic  moduli 

Isotropic  moduli  defined  in  eqn.  (B-5) 

Second  and  third  order  moduli  in  Brugger's  notation 
Isotropic  moduli  in  Voight's  notation 
Elastic  wave  speeds  in  isotropic  materials 
Defined  in  eqn.  (4.13) 

Unit  vectors  along  X,  Y  axes 

Deformation  gradient 

Deformation  gradients  at  the  boundary 

Unit  proper  vector  in  simple  wave 

Unit  normal  vector 

Order  of  magnitude  symbols 

Unit  proper  vectors  for  linear  wave 

Acoustic  tensor 

Stress  tensor 

Unit  tangent  vector 

Time 

Particle  velocity 
Wave  speed 

Linear  elastic  wave  speeds,  defined  as  in  (2.6) 

Steady  speed  parallel  to  the  boundary 


IX 


V 


o,  ao,  a,  a 


5,  n,  c 


♦  ,  <!»>  to 


Perturbation  of  shock  speed 
Strain  energy  function 
Material  coordinate 
Spatial  coordinate 
Deformation  gradient 
Small  angles 
Simple  wave  parameter 
Variation  across  i**1  wave 
Strain  tensor 
Wave  angles 
Limiting  angle 
Lame  constants 
Components  of  a 
Density 

Second  order  corrections  to  amplitudes 
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Prof.  A.  Pipkin 
Prof.  P.  Symonds 
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East  Lansing,  MI  48823 

1  New  York  University 

Department  of  Mathematics 
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Dr.  W.  Baker 
8500  Culebra  Road 
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ATTN:  Prof.  R.M.  Haythomthwaite 
Philadelphia,  PA  19122 

1  Tulane  University 

Dept  of  Mechanical  Engineering 
ATTN:  Dr.  S.  Cowin 
New  Orleans,  LA  70112 

2  University  of  California 
ATTN:  Dr.  M.  Carroll 
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USER  EVALUATION  OF  REPORT 


Please  take  a  few  minutes  to  answer  the  questions  below;  tear  out 
this  sheet,  fold  as  indicated,  staple  or  tape  closed,  and  place 
in  the  mail.  Your  comments  will  provide  us  with  information  for 
improving  future  reports. 

1 .  BRL  Report  Number _ 

2.  Does  this  report  satisfy  a  need?  (Comment  on  purpose,  related 
project,  or  other  area  of  interest  for  which  report  will  be  used.) 


3.  How,  specifically,  is  the  report  being  used?  (Information 
source,  design  data  or  procedure,  management  procedure,  source  of 
ideas,  etc.) _ _  _ _ 


4.  Has  the  information  in  this  report  led  to  any  quantitative 
savings  as  far  as  man-hours/contract  dollars  saved,  operating  costs 
avoided,  efficiencies  achieved,  etc.?  If  so,  please  elaborate. 


5.  General  Comments  (Indicate  what  you  think  should  be  changed  to 
make  this  report  and  future  reports  of  this  type  more  responsive 
to  your  needs,  more  usable,  improve  readability,  etc.) _ 


6.  If  you  would  like  to  be  contacted  by  the  personnel  who  prepared 
this  report  to  raise  specific  questions  or  discuss  the  topic, 
please  fill  in  the  following  information. 

Name : 


Telephone  Number ; 
Organization  Address; 
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US  Army  Ballistic  Research  Laboratory 
Aberdeen  Proving  Ground,  MD  21005 


OFFICIAL  BUSINESS 

PENALTY  FOR  PRIVATE  USE.  S300  BUSINESS  REPLY  MAIL 

FIRST  CLASS  PERMIT  NO  12062  WASHINGTON, DC 
POSTAGE  WILL  BE  PAIO  BY  DEPARTMENT  OF  THE  ARMY 


Director 

US  Army  Ballistic  Research  Laboratory 
ATTN :  DRDAR-TSB 

Aberdeen  Proving  Ground,  MD  21005 
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